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The suspension noise in interferometric gravitational wave detectors is caused by losses at the top 
and the bottom attachments of each suspension fiber. We use the Fluctuation-Dissipation theorem 
to argue that by careful positioning of the laser beam spot on the mirror face it is possible to 
reduce the contribution of the bottom attachment point to the suspension noise by several orders of 
magnitude. For example, for the initial and enhanced LIGO design parameters (i.e. mirror masses 
and sizes, and suspension fibers' lengths and diameters) we predict a reduction of ~ 100 in the 
"bottom" spectral density throughout the band 35 — 100Hz of serious thermal noise. 

We then propose a readout scheme which suppresses the suspension noise contribution of the top 
attachment point. The idea is to monitor an averaged horizontal displacement of the fiber of length 
I; this allows one to record the contribution of the top attachment point to the suspension noise, 
and later subtract it from the int erferomet er readout. This method will allow a suppression factor 
in spectral density of 7.4 (l/d 2 ) yj Mg/nE, where d is the fiber's diameter, E is it's Young modulus 
^— i , and M is the mass of the mirror. For the test mass parameters of the initial and enhanced LIGO 

designs this reduction factor is 132 x (Z/30cm)(0.6mm/d) 2 . 

We offer what we think might become a practical implementation of such a readout scheme. We 
propose to position a thin optical waveguide close to a fused silica fiber used as the suspension fiber. 
^ ' The waveguide itself is at the surface of a solid fused silica slab which is attached rigidly to the 

last mass of the seismic isolation stack (see Fig. 5). The thermal motion of the suspension fiber is 
recorded through the phaseshift of an optical wave passed through the waveguide. A laser power of 
lmW should be sufficient to achieve the desired sensitivity. 
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CT 1, I. INTRODUCTION 
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Random thermal motion will be the dominant noise source in the frequency band of 35 — 100 Hz for the first inter- 
ferometers |jj and in the frequency band of 25 — 126 Hz for the enhanced interferometers [] in the Laser Interferometer 
Gravitational Wave Observatory (LIGO) 0. 

The thermal noise in this frequency band is caused by the losses in the suspension fibers, in particular at the top 
and the bottom of each fiber's attachment point. So far the only known way to reduce the thermal noise has been to 
improve the quality of the suspension fibers and their attachments. Here we suggest a different approach: 

In Section II we will present a general analysis of the suspension noise based on a direct application of the 
Fluctuation-Dissipation theorem. We will explicitly separate the contributions to the thermal noise of the top and 
the bottom attachment points of the suspension fibers. It has been a common opinion that the top and bottom 
attachments contribute equally to the thermal noise. We shall challenge this point of view. In fact, we will show 
that if one shifts the laser beam spot down from the center of the mirror by an appropriately chosen distance h, the 
contribution of the bottom attachment point to the thermal noise can be reduced by several orders of magnitude. 
Fig. 3 presents plots of this reduction factor in the frequency band 35-100Hz for three different choices of h. What is 
plotted here is the ratio 5bottom(/)/<Stop(/), where «Sbottom(/) and <Stop(/) are the spectral densities of thermal noise 
contributed by the bottom and the top attachment points respectively. All three values of h are close to 



1 To be specific, we refer to the step 4 of LIGO enhancement — see In these the suspension thermal noise was calculated 
assuming the structural damping mechanism. However, the nature of dissipation in fused silica (e.g. viscous vs structural) is 
not yet fully established for the above frequency bands. 

2 The analysis of this paper is fully applicable to all other Interferometric Gravitational Wave detectors (e.g. VIRGO, GEO-600, 
TAMA etc.). For the sake of brevity in this paper we will refer only to LIGO. 
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[cf. Eq(|lJ)], where I is the length of the suspension fiber, / is the test- mass moment of inertia for rotation about the 
center of mass in the plane of Fig. 1 (see later), R is the radius of the mirror face and M is the mass of the test mass. 
The numerical values of these parameters for the initial and enhanced LIGO interferometers are 

M = 10kg, I = 30cm, R = 12.5cm, (2) 



I = 4.73 x 10 5 g cm 2 , h = 1.11cm. 

Out of the three graphs presented in Fig. 3 the one with h = 1.0cm seems to be the optimal one. From the graphs 
we see that reduction factors of ~ 10 -2 in the "bottom" component of the thermal noise is possible over the entire 
band of serious thermal noise: 35 to 100 Hz. 

In Sec. IIIA we concentrate on the top attachment point. Lossy defects at the top create noise not only in 
the test mass motion, but also noise in the motion of the fiber. The latter is significantly larger than the former 
- by a factor of order / 2 //p en d u ium a ^ frequencies above the pendulum frequency and below the violin resonances 
(which are the frequencies of interest for LIGO thermal noise). We show that if one monitors the average horizontal 
displacement of the suspension fiber of length I, one can essentially record the fluctuating "driving force" originating at 
the suspension top, and then subtract it from the interferometer's readout, thereby reducing thermal noise originating 
at the suspension top. The reduction factor in the spectral density of thermal noise is given by P — 0.93 ■ l/X [cf. Eq( 
!§)]. Here 

A = (d 2 /8) y / irE/Mg (3) 

is the length of the segment of fiber near it's top where the bending is greatest, d is the fiber's diameter, E is the 
fiber's Young modulus and g is the acceleration of gravity. For a fused silica fiber of diameter d = 0.6mm one gets a 
thermal noise reduction factor of P ~ 132. 

In Sec. IIIB we offer a particular way of implementing such a procedure. The basic idea is shown in Fig. 5. A fused 
silica slab is rigidly attached to the "ceiling" (i.e. to the last mass of the seismic isolation stack), and a waveguide ab 
is carved into the slab's surface. A monochromatic optical wave is set up in the waveguide, and a fused silica fiber 
used as the suspension fiber is positioned close to the waveguide, within the optical wave's evanescent field. When 
the fiber is displaced relative to the waveguide, it will change the optical wave's propagation speed, thus inducing an 
overall phaseshift of the wave. The detailed calculations in Sec. IIIB show that ~ lmW of optical power in the wave 
is sufficient to reach the required sensitivity. 



II. HOW TO REDUCE THERMAL NOISE ORIGINATING AT THE BOTTOM ATTACHMENT POINT 



A. The model and formalism 



The particular suspension that we consider is sketched in Fig. 1. We consider a compact rigid test mass of mass M 
suspended by a single fiber of length I and mass m; the fiber's bottom end is attached, for concreteness, to the top of 
the test mass (the main conclusions of this paper are also valid when the test mass is suspended by a fiber loop, as is 
planned for LIGOV. 

References j|, Q|, f|, give detailed explanations of how to use the Fluctuation-Dissipation theorem directly 
(without normal-mode decomposition) to calculate the spectral density of thermal noise ^. In what follows we use 
the approach elaborated in . 

To calculate the spectral density S K (f) of suspension's thermal noise at frequency / we imagine applying an 
oscillating force F perpendicular to the test mass's mirror surface at the center of the readout laser beam spot Ft 



3 The original formulation of the Fluctuation-Dissipation theorem is given in |j 

4 This prescription is only valid when the test masses are perfectly rigid, which is a good approximation when dealing with 
suspension thermal noise. The case when the test masses are no longer considered to be rigid (e.g. for an internal thermal 
noise calculations) is treated in detail in In that case the force F(t) must be spread out over the laser beam spot instead 
of applied to it's center point. 
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F(t) = F cos(2irft). (4) 

Then S x is given by [cf. Eq (3) of §] 

2k B T W di: 
* 2 P F$ 



S *(f)= J f2 T , (5) 



where Wdi ss is the average power dissipated in the system (suspension, in our case) when the force F(t) is applied, 
fee is Boltzmann's constant and T is the temperature. 

For concreteness, assume that the dissipation in the fiber occurs through structural damping (our conclusions 
will hold equally well for viscous or thermoelastic damping). In this case, the average power dissipated during the 
oscillatory motion of frequency / is given by Q 

W diss = 27r/C/ max <H/), (6) 

where J7 max is the energy of the fiber's elastic deformation at a moment when it is maximally bent under the action 
of the oscillatory force in Eq. (|j), and 4>{f) is the "loss angle" of the material. The energy of the fiber's elastic 
deformation is given by 



U= J -^j\z[y»] 



(7) 



where E is the Young modulus of the fiber material, J is the geometric moment of inertia of the fiber (for a fiber with 
circular cross section of diameter d one has J = ird 4 /64), z is distance along the fiber with z — at the top and z — I 
at the bottom, and y(z) is the fiber's horizontal displacement from a vertical line. 

This method of calculating thermal noise is useful for a qualitative analysis of the system, as well as quantitative 
analysis. In particular, it allows one to see which part of the suspension fiber contributes the most to the thermal 
noise. Assume, for a start, that the laser beam is positioned exactly in the middle of the mirror. Then to work out the 
thermal noise one has to imagine applying the oscillating force F in Eq. (3) to the mirror center; the motion of the 
fiber and the mirror under the action of the force are shown in Fig. 2a. Here we assume that the detection frequency 
/ (and hence the frequency of the applied force) satisfies f p , f r << f << f v , where f p , f r , f v are the frequencies of 
the pendulum, rocking and first violin mode respectively (this condition implies that horizontal and rotational motion 
of the test mass is not affected by the presence of the fiber, and that the fiber itself remains straight). 

From Fig. 2a it is clear that the fiber bends equally at the top and the bottom (we always assume that at the 
attachment point the fiber has to be normal to the surface to which it is attached). The total energy of elastic 
deformation is 

2 _ MgA f F 



where A = y JE/Mg is the characteristic length over which the fiber is bent near the attachment points, uj = 2irf is 
the angular frequency of detection, and a is the angle between the straight part of the fiber and the vertical. 

The bending of the fiber at the bottom can be avoided if one applies the force F in Eq. (^) not at the middle of the 
mirror, but at some distance h below the center. In particular, we should choose h so that the mirror itself rotates by 
the same angle as the fiber under the action of the applied force; the resulting motion is shown on Fig. 2b. Physically 
this means that if we position our laser beam at a distance h below the mirror center, then the bottom attachment 
point will not contribute to the thermal noise when h is carefully chosen. This means that the overall suspension 
noise will be reduced by a factor of order 2 (in fact, more precisely, by a factor of 2(1 + R/l), where R is the radius 
of the mirror and I is the length of the string, — see later in this section) . 

In the rest of this section and Appendix A we find the general expression for the suspension thermal noise, and 
we then work out the optimal h for the frequency band of interest for LIGO. We will assume that when a periodic 
oscillation of frequency / is induced in the system, the average power dissipated as heat in the suspension is given by 

Wd iss = / [Ctop(/)4 + Cbottom (/)<*!] • (9) 

Here 5t and <3b are the amplitudes of oscillations of the angles ax and «b respectively (see Fig. 1), and Ctop an d 
Cbottom are frequency-dependent quantities characterizing dissipation at the top and the bottom respectively. For the 
case of structural damping 

Ctop = Cbottom = nf<f>(f)Mg\, (10) 
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where A is given by Eq. ([|) of the introduction. 

To compute Wdiss we need to evaluate <St and <5b by analyzing the dynamics of the oscillations. This is done in 
Appendix A, see Eqs. ( flO| ) and (|39|). Putting these equations into Eq. (g) and then into Eq. (JsJ) , we obtain [cf. Eq. 
©] 



S*(J) 



Ak B T f I/M - R(g/u 2 + h) ' 

ttcj \ [i ff - MgR {g/uj 2 - R)} cos(fcZ) - (iw 2 - Af pi?) sin(fcZ) /fc 



Ctop + Cbottom COS 2 (kl) 



I/M -h[R + tan (kl) /k - g/u 2 } 
I/M - R(g/to 2 + h) 



(11) 



Here k = u>/c — 2nf/c, c — y/glM/m is the speed of propagation of a transverse wave in the fiber. From the above 
equation we can infer the ratio of the bottom and the top contributions to the thermal noise: 



iSbottom(/) Cbottom(/) 



StoM) 



Ctop(/) 



cos 2 (kl) 



I/M -h[R + tan (kl) /k - g/u 2 
I/M - R(g/to 2 + h) 



(12) 



This is the most important equation in this section of the paper; it will be discussed in the next subsection. 



B. The case of low- frequency suspension noise 



When the detection frequency / is far below the frequency of the fundamental violin mode, f v , then kl <C 1 in Eq. 
(U) and 



tan(kl) 



I 



l + \{klf 



(13) 



Let us assume that the top and the bottom are equally lossy, i.e. Ctop = Cbottom, as the would be for structural 
damping, Eq (|io| ) above. We choose h to be 



I 



M(R + 1) ' 



Putting Eqs @ and @ into Eq. @, we get 

<Sbottom(/) 7T 4 



Sto P (/) 9 [i_ (R/h)(uj/u 2 )] 2 V/v 



(14) 



(15) 



where uj p — \fgfl. 

For the initial and enhanced LIGO design f v ~ 400Hz, M ~ 10kg, I ~ 4.73 x 10 _2 kg x m 2 , R ~ 12.5cm, and the 
interesting frequency range where suspension noise is expected to dominate is 35 — 100Hz (actually, this depends on 
the stage of enhancement. The frequency band specified above is where the suspension thermal noise is expected to 
dominate in the initial LIGO; in the enhanced version this frequency interval will be larger). In this case Eq. ( |l5| ) 

gives Sbottom(/)/£top(/) " °- 002 ~ °- 2 - 

In Fig. 3 we give plots for <Sbottom/Stop as a function of the detection frequency / for three different choices of h. 
We have used Eq. (|l5|) to make all the plots and we set i, M , R and I to the numerical values appropriate for the 
initial and enhanced LIGO design and given at the beginning of this section. 

The first curve is plotted for h given by Eq. (|i~4j), in our case h = 1.11cm. The second and third curves are for 
h = 1.0cm and h — 0.9cm ; these values of h are chosen so that S'bottom/S'top = for / = 80Hz and / = 105Hz 
respectively. Out of the three cases the choice h = 1cm gives the best overall performance across the considered 
frequency band, with the typical reduction factor of 



She 



10- 



top 



(16) 



^From Eq. (|TT]) we see that choosing h close to the value in Eq. ( |l4| ) reduces the total suspension thermal noise by a 
factor close to 2(1 + R/l) ~ 3 relative to the case when h = 0. 
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C. High-frequency suspension thermal noise 



A somewhat less interesting observation is that for h = and / n = f v (n + 1/2), where n is an integer, 

5'bottom(/n) 



Stop(/n) 



0. (17) 



Unfortunately, at / = f n the interferometer's noise is dominated by shot noise. However, if one uses an advanced 
optical topology — for example, resonant sideband extraction — then it is possible to reduce the shot noise in a 
narrow band around any chosen frequency. Then the thermal noise may dominate in this narrow band, and our 
observation ( p~7|) may be useful in case one tries to reduce the thermal noise by cooling of the fiber top. 

III. HOW TO CONTROL NOISE FROM THE TOP 
A. The concept 

In this section we propose a recipe for how to decrease the influence of the thermally fluctuating stress at the top 
part of the suspension fiber. The basic idea is the following: 

Intuitively, the fluctuations at the top cause bending of the fiber at the top, which will be a random process in 
time. This random bending will randomly move the rest of the fiber and ultimately drive the random motion of the 
test mass. We propose to measure directly the thermally driven fluctuations in the horizontal displacement of the 
fiber, and from them infer the fluctuating force which drives the random motion of the mirror. We can then subtract 
the motion due to this fluctuating force from the interferometer output^]. 

Formally this amounts to introducing a new readout variable q as follows: 

(} = -^mirror + -^fibei- (18) 

Here X m ; rror is the horizontal displacement of the laser spot's center (i.e. the signal ultimately read by the interfer- 
ometer's photodiode), and 

X nhel = / dz<£(z)y(z) (19) 
Jo 

is the fiber's horizontal displacement weighted by some function &(z) to be discussed below. We will postpone the 
discussion of how to measure q experimentally until the next section; here we concentrate on finding the optimal $(z) 
and seeing what is the maximal possible reduction in the thermal noise. 

To find the spectral density of fluctuations in q we need to imagine acting on the system with sinusoidal force 
F q oc cos(2-7r ft) that appears in the interaction hamiltonian in the following way 



Hint = 



-qF q = -X milIOI F q - [ dzF q $(z)y(z); (20) 
Jo 



cf. the discussion of the Fluctuation-Dissipation theorem in Ref . . From the Eq. ( pCf ) we observe that applying the 
generalized force F q to the system is equivalent to applying two forces simultaneously: one is a force of magnitude F q 
applied to the mirror surface at the center of the beam spot, and the other is a force distributed along the fiber in 
the following manner: 

^ = F 9 *(z). (21) 

The resulting motion of the system is shown in Fig. 4. The intuitive idea is to choose the weighting function $>(z) so 
that when the beam spot's height h has also been appropriately chosen, F q induces no bending of the fiber at the top 
or at the bottom. 



The idea of thermal noise compensation is not new (e.g. @], [0). However, our detailed treatment and concrete experimental 
proposal is different from anything prior to this paper. 



5 



In the case of structural damping the dissipated power is proportional to the elastic energy U of the fiber. Thus 
formally one has to choose $>(z) and h so that U is minimized. It is convenient to reformulate the problem: to find 
the shape of the fiber y(z) and beam-spot height h for which the functional in Eq. (^) has a minimum, and after this 
calculate the distribution $(z) of the driving force on the fiber that will produce the desired shape y(z). In Appendix 
B we carry out this straightforward but somewhat tedious task. We obtain [cf. Eq. (|5C 



F / z \ * 

2/optimal(20 = (j) 



F„ 



Mlu 2 



3{r + l)-z/l 
2(3r 2 + 3r + l) 

(t)> 76 -°' 18 t 



(22) 



Here r = R/l, R is the radius of the mirror, I is the length of the fiber, u> = 2irf is the angular frequency of detection. 
We substitute here and below r = 0.42 corresponding to the initial and enhanced LIGO test masses. The profile of 
the distributed force acting on the fiber and hence of <&(z) is mainly determined by y"(z) (see Appendix B): 



$(z) ~ $ 



M 9 n, \ 

-f^y (*)> 



(23) 



which gives [cf. Eq. (53)] 



*o(*) - 

u 2 l 



z 



3r 2 + 3r + 1 



1.53 



1.08? 



(24) 



where ui p = \Jg/l. When the force distribution has this optimal form, the elastic energy has the minimum value 



XMgX ( F q 



3r 2 + 3r+ll 2 \Mw 2 l 

1.08A TT 

x U , 



I 



(25) 



where Uq is the elastic energy in Eq. (||). Therefore, for a fused silica fiber with E ~ 6.9 x 10 10 Pa and d = 0.6mm, 
we get A ~ 2.1mm and the maximal reduction factor for the spectral density of suspension thermal noise is 



P = 



I 

1.08A 



132. 



(26) 



B. Experimental realization— a proposal 

1. Preliminary remarks 

Before describing a particular experimental realization of the above scheme, a few general remarks are in order. 

First, one might worry that our averaging function $(z) is frequency dependent — in general, that could make the 
experimental implementation very difficult. In particular [see Appendix B, Eq.([52|)], $ consists of two components: 
$ = $o + $1, where $o an d $i as given by Eq. (|5^) have very different frequency dependence. However at the 
frequencies of interest <i>o ^> $i, and then the approximate formula ( |24| ) for the averaging function $(z) = <&o{z) 
is a product of two terms: one which depends only on the frequency / (i.e. $(2) oc l// 2 ), and the other which 
depends only on the coordinate z. This feature makes the scheme feasible for a broad range of frequencies. It is 
sufficient that our device measures the displacement of the fiber with the frequency-independent averaging function 
$(z) oc f 2 x 3>(z), and that the frequency dependence is then put back in during data analysis when constructing the 
readout variable q: 

q = A mirror + T)(f) [ dz$(z)y{z), (27) 
Jo 

where r](f) oc /~ 2 is chosen so that 77$ = $. 







As mentioned above, Eq. ( p4| ) is an approximation valid when the fiber has no inertia, i.e. when / <C fv = 
(lowest violin-mode frequency). When the inertia of the fiber becomes important (<J>i ~ $o)j it is no longer possible 
to factor out a frequency-dependent part of As a result, when / gets closer to /„, the effectiveness of the thermal 
noise suppression (i.e. the value of P) is reduced. A detailed analysis shows that if we choose $(z) so that the thermal 
noise compensation is optimal (P = P max ) at low frequencies / -C f v , then at / = 0.2/ u we have P ~ 0.9P max , at 
/ = 0.32/u we have P ~ 0.5P max , and beyond this P is reduced sharply as we approach the first violin mode. For 
the fused silica fiber discussed above f v ~ 400Hz, so the compensation is effective throughout the band 35 — 100Hz 
where suspension thermal noise dominates. It is worth emphasizing that this deterioration in the reduction factor 
only happens when we use the averaging function <&o instead of $0 + $1 close to the violin frequency. Thus, this 
limitation is one of technology and not of principle. Perhaps, it is possible to conceive of a scheme where the correct 
averaging function is implemented at all frequencies. However, we have not been able to do so. 

Secondly, any sensor used for monitoring the fiber coordinate A^bcr wn l have an intrinsic noise which will deteriorate 
the quality of the thermal-noise compensation. In particular, the overall reduction factor P c g is given by 

1 _ _j_ _|_ *Sfiber meas(/) (28) 
PcS P S R bcr therm 

(/)' v ' 

where S^ber meas(/) is the spectral density of intrinsic noise of the device which measures the average displacement 
of the fiber and Ssb el therm(/) is the spectral density of thermal fluctuations of the same displacement. 
For the case of structural damping it is easy to estimate 



VS X fiber thorm(/)/ ~ J ^ ~ lO^Cm, (29) 

where we assume that <f> ~ 10 -7 for fused silica. If our goal is to achieve P ~ 100 then the condition P g- ~ P implies 



V^fibcr mcas/ « Sj 5fibG1 ^ Crm/ ~ 10- 15 Cm. (30) 

We shall take the above number as a sensitivity goal that our measuring device should achieve. 



2. Proposed measuring device 

Now we are ready to describe a possible practical implementation of our thermal-noise compensation scheme. Figure 
5 illustrates the basic idea. We propose to use a fused silica optical fiber with the refractive index n\ for the test 
mass's suspension. Next to this fiber we attach to the top seismic isolation plate (i.e. the "ceiling") a rigid block 
of the fused silica A with the same index of refraction ri\. On the surface of this rigid block we put a thin optical 
waveguide with refractive index n% such that n-i > n\, so that the waveguide is at a distance ~ A op ti C ai/27r from the 
suspension fiber. It is assumed that the side of the waveguide close to the suspension fiber does not have any coating, 
i.e. it is "naked" . In this configuration the optical wave may propagate through the waveguide without substantial 
scattering even though the suspension fiber is within the wave's evanescent zone. This device will produce a relatively 
large response to the displacement Afiber in the form of a phaseshift of A<fi of the optical wave: 

A0 = X 27rXfibCT 2nl , (31) 

^optical ^optical 

where the dimensionless factor K depends on the values of n\ and ri2 and for typical optical waveguides is K ~ 10 -3 . 
Equation (|3l|) implies that in order to register Xsb el ~ 10~ 15 cm we need a sensitivity A(f> ~ 10 -7 . Thus for averaging 
time of Tg ra v = O.Olsec we need to use the power of coherent light of W ~ lmW. This power can be decreased if one 
uses a resonant standing wave in the waveguide. 

Apart from the shot noise of the laser light, let us briefly discuss two other kinds of noise in this sensor. A more 
complete discussion will be presented elsewhere. 

The first kind is seismic noise. A simple calculation shows that the seismic contribution to the noise in the readout 
variable q is about twice as large in spectral density as the seismic contribution to the noise in A m i rror . Thus the 
seismic noise will not be an issue at frequencies above the "seismic wall" of the LIGO sensitivity curve. 

The second kind of noise we want to mention is the mechanical thermal fluctuations of the waveguide itself. Our 
estimates show that if these fluctuations are caused by structural damping (and not by some surface or contact 
defects), then the ratio of the mechanical thermal fluctuations of the waveguide to those of the fiber is 
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Thus, if the system is sufficiently clean then the mechanical thermal fluctuations of the waveguide will probably not 
significantly reduce the sensitivity of our sensor. 

It is worth noting that in order to achieve the optimal compensation of thermal noise, the distance d(z) between 
the suspension fiber and the waveguide has to vary in accord with the optimal profile of the averaging function: 

d = A-B\og[${z)], (33) 

where A and B are constants to be discussed elsewhere. In this case the phase of the waveguide's output records the 
optimally averaged coordinate Xgber of the fiber. 

The profile d(z) may be difficult for experimental realization. However we find that in the simplest case when $(z) 
is a constant over the length I of averaging, the factor P is reduced very little: from P = 132 to P ~ 120. 



IV. CONCLUSION 



In this paper we have done two things. 

Firstly, we have shown that by an appropriate positioning of the laser's beam spot on the surface of each test-mass 
mirror, one can reduce the contribution of the suspension fiber's bottom to the suspension thermal noise by two to 
three orders of magnitude in the frequency band of 35 — 100Hz for the initial LIGO design. 

Secondly, we have proposed a way to compensate the suspension thermal noise originating from the top of each 
fiber by monitoring independently the fiber's random horizontal displacement. In the best case, with the system 
parameters for the initial or enhanced LIGO design, one can get a reduction factor of the order of P = 130 in spectral 
density over the entire 35 — 100Hz band, when both the first and second procedures are applied; and with realistic 
defects in the design one should be able to get a reduction of at least P ~ 100 

The device that compensates the suspension thermal noise can ease the requirements to quality of suspension 
system. In particular, if this device allows the reduction factor of P — 100, this would effectively increase the 
quality factors of pendulum and violin modes by a factor of P = 100. So far the highest quality factor Q ~ 10 s of 
the pendulum mode was achieved in Q for a fused silica suspension fiber, which allows one to reach the Standard 
Quantum Limit for averaging time of 10 _3 sec. Implementation of our proposal could effectively increase this quality 
factor to Qcff ~ 10 10 , which would reduce the thermal noise in LIGO to the level of Standard Quantum Limit for 
averaging time of 10 _2 sec. Then the techniques which allow one to beat the Standard Quantum Limit (see e.g. p3|) 
could be used in the enhanced LIGO interferometers. 
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APPENDIX A 



In this appendix we solve the dynamical problem of finding the amplitudes c*t and cub of oscillation of the top and 
bottom bending angles in Eq. (^|) when a periodic force 

F = F cos(ujt) (34) 

is applied to the mirror at a distance h below the mirror center [we use these amplitudes in Eq. of the text]. For 
convenience we complexify all of the quantities: 

F77T „lU)t „ — „lLUt „ — „liilt 

= roe , qt = CKTe , cub = c^Be , 



tp = 4>e l 



where x is the displacement of the test mass's center of mass and ip is the angle by which the mirror is rotated (see 
Fig. 1) under the action of the force F(t). As usual, ui = 2nf is the angular frequency. 
^From the projection of the Newton's Second Law on the horizontal axis we have 



F - (a B - ip)Mg = -Mu 2 x, 
and, for the rotational degree of freedom, the equation of motion is 

F h + MgRotB = Iuj 2 ^, 



(35) 



(36) 



where R is the radius of the test-mass cylinder and I is the moment of inertia for rotation about the test-mass center 
of mass in the plane of the Fig. 1. In the two equations above we assume that «b and tp are small. 
The fiber's horizontal displacement y from a vertical line approximately satisfies the wave equation: 



d 2 y 
8t 2 



" dz 2 ' 



(37) 



where z is distance along the wire, with z = at the top and z = I at the bottom, and c = y/glM/m is the transverse 
speed of sound in the wire. In this Appendix we use Eq. fl37j ) for flexible wire since it's solutions are simple. If one 
takes the stiffness into account this changes the solutions of Eq. (|37j ) by a relative order of A/7, see e.g. jjjl. However, 
when using Eq. (|37|), we must allow non-zero bending angles at the top and bottom attachment points, ax and a B - 
The energy of elastic strain of the wire then consists of two components: one from the bulk of the wire given by Eq. 
(0) , and the other from the bending at the attachment points given by Eq. (||) . The solution to Eq. (||^) is 



y(z,t) = A sin (fez) e 4 ' 



(38) 



where k = uj/c is the wave vector of an off-resonance standing wave induced in the fiber and A is a constant. The 
boundary condition is set at the bottom by 



A sin (kl) 
fcAcos (kl) 



x + Rtjj 



Putting these two equations into Eqs. ([ 

-F 



and 



qt = —Fq-. 



5J) and (§|), we find 

I/M - h [R + tan (kl) /k - g/u 2 ] 
MgR 2 + (Iuj 2 - MgR) [g/uj 2 - tan (kl) Jk] 

I/M - R(g/uj 2 + h) 



(39) 



[Ig - MgR (g/cu 2 - R)] cos(fcZ) - (Iuj 2 - MgR) sin(kl)/k 
Putting Eqs. ([59]), (40) and (|9|) into Eq. (||), we finally get for the spectral density of the suspension thermal noise 



(40) 



s,(f) 



4k B T 



I/M - R(g/u 2 + h) 



iruj { [Ig - MgR (g/uj 2 - R)] cos(kl) - (Iuj 2 - MgR) sin(kl)/k 



Ctop + (bottom COS 2 (kl) 



I/M -h[R+ tan (kl) /k - g/uj 2 
I/M - R(g/u 2 + h) 



-i 2 



(41) 



APPENDIX B 



Here we calculate the optimal shape y ptimai(.z) of the fiber and the vertical position of the laser beam spot h that 
minimize the fiber's elastic deformation energy [Eq. (Q)]. 

It is easy to deduce from Eq. (^) that energy minimizing function y(z) obeys the equation y""(z) — 0. Therefore 

y(z) z z 2 % 

'—j—=a + a 1 j + a 2 j I + a 3 -jj, (42) 
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where are constants to be determined. 

Let us discuss the boundary conditions. Strictly speaking, the boundary conditions should be such that the 
fiber is perpendicular to the surface of attachment at both the top and the bottom. Therefore at the top we have 
y(0) = y'(0) = 0, from which immediately follows ao = a\ = 0. However, at the bottom it is more convenient for 
our calculations embody the bending of the fiber, on the lengthscale A, in a bending angle a B as in Fig. 1, and 
correspondingly add an additional term 

t/add = (l/4)Af 5 Aa| (43) 

to the energy functional in Eq. (Q), and then in Eq. (^2|) evaluate y(l) and it's derivatives above the A-scale bend. 
Our energy minimization procedure will make the angle a B so small that the additional elastic energy as given by 
Eq. (ff3| ) is negligible compared to U in Eq. (0) 

The coefficients ai and 0,3 can be inferred from force and torque balance at the test mass: 

F q - Mgy'{l) = -Mw 2 (y(l) + R(y'(l) + a B )), (44) 

and 

F q h - MgRa B = -Iuj 2 (y'(l) + a B ). 

It is useful to rewrite these equations in a dimcnsionless form: 

£(1 + r)(r - a)) + ra B = -£o, 

v£ + a B(l — M rc = —fis£ ; (45) 



where 



yJJ) _ y'(l)l _h 

* 1 ' 11 y(i) ' 3 r 



K , „ Ml 2 

a = -§ ~ 1CT 3 1Q- 6 , r = 0.42, fi = —— = 19, 



where to p = yg/l- Here we have used for estimates the mirror parameters for the initial and enhanced LIGO 
interferometers. Solving the above system of equations (Mq) for £ and a B (taking 77 as a parameter) we get: 



„ 1] — LbsCl + T](r — a)) 1 \m 

"B = £07— — Try- ; ~ £0 7? - HS{1 + r) r - a 

[1 + 7](r — a)\[l — fxra] — rq 

£ = -go h - , vim i ~ V r ( a + s )) 

[1 + rj[r — a)\[l — fira\ — rr\ 

Let us choose the parameter s so that a B — for some angular frequency ujq in the frequency band 35 — 100Hz where 
thermal noise is most serious: 



fi[l + rj(r - a )} + r/r] ' uj 2 



, a = ^ (46) 



Then we get for a B and £ 



v 2 

a B ^£,o— (a-ao) (47) 

1 + rjr 

1 



1 + ijr 



We can express the coefficients 03 and a 2 in terms of £ and rj by combining Eqs. ( |4^ ) and (|45|), and we can then 
calculate the elastic energy according to Eq. (0) : 
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\Muj 2 lJ I (1+rr?) 



This function has the minimal value 



at optimal r\ given by 



I 3 1.08A 



Here C/o is the energy of elastic strain of the fiber when the force of magnitude F q is applied in mirror center, as 
worked out in Eq. (ph. Now we can figure out the optimal shape of the fiber's horizontal displacement: 

F q fz\ 2 ( 3(r + l)-2// 



2/optimal(2) - ( l ) 



2(3r 2 + 3r + l) 



F, 



0.76-0.18- . (50) 



Muj 2 \IJ V I 

From Eq. (|4q ) we get h = I x s ~ 1.55cm. 

Using (|4T|) one can show that «b < 1-7 ■ 10~ 3 • £0 over the frequency band 35 — 100Hz. From this and Eq. fl43|), one 
can compute the energy due to the bending at the fiber bottom: [7 a( jd — 1-4 • 10~ 6 x Eq. We see that [/ a dd <C C/ m i rl 
and hence over the frequency band of interest the small bending at the bottom does not contribute significantly to 
the total energy of elastic deformation. 

The profile of the distributed force and correspondingly the function $ are given by 

- - P io 2 y(z) - Mgy"{z) + IEy""(z). (51) 

Here p is the fiber density per unit length. Since y""{z) — 0, the function $ consists of two terms Q(z) — $ Q (z) + $>i(z), 
where 

$ (z) = -^V(*), 9l {z) = -^-y{z). (52) 



$oW _ ± . (! + r _ |) . . . ( L53 _ !.„«£) . (53) 

We see that $0 is much greater than $1 in our frequency range (10 — 100Hz for the initial LIGO). 
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FIG. 1. We consider a test mass suspended on a single fiber. The fiber's bottom is attached to the top of the test mass, and 
the fiber's top is attached to the last stage of the seismic isolation stack. It is assumed that at attachment points the fiber is 
perpendicular to the surface to which it is attached. 





FIG. 2. Motion of the test mass and the suspension fiber under the action of an oscillating force applied at the center of the 
laser beam spot in two different cases: a) the beam spot is positioned at the mirror center, the fiber bends equally at the top 
and the bottom, and b) the position of the beam spot is shifted down from the center of the mirror, so that there is no bending 
of the fiber at the bottom. 
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FIG. 3. Fig. 3. A plot of Sbottom(f) / St op (f) as a function of frequency / for three different positions of the laser beam spot. 




FIG. 4. Motion of the test mass and suspension fiber under the action of the generalized force F q defined in Eq. ( J20| ) of the 
text. The force F q should be chosen so that there is no bending of the fiber at it's top and bottom attachment points. 




waveguide 



suspension wire 

FIG. 5. A proposed scheme for compensation of a suspension thermal noise. The optical waveguide ab is positioned close to 
the suspension fiber made of fused silica. A horizontal displacement of the suspension fiber is recorded through a phase shift 
of an optical wave propagating through the waveguide. 
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